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Abstract 

The null horizon focusing equation is equivalent via the fluid/gravity correspondence to the 
entropy balance law of the fluid. Using this equation we derive a simple novel formula for the 
bulk viscosity of the fluid. The formula is expressed in terms of the dependence of scalar fields 
at the horizon on thermodynamic variables such as the entropy and charge densities. We apply 
the formula to three classes of gauge theory plasmas: non-conformal branes, perturbations of the 
A/" = 4 super symmetric Yang-Mills theory and holographic models of QCD, and discuss its range 
of applicability. 
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I. INTRODUCTION AND SUMMARY 

One important application of the AdS/CFT correspondence is the calculation of 

transport properties such as viscosity and conductivity of strongly coupled large Nc gauge 
theory plasmas using methods from black hole physics. This follows from the map of thermal 
states of field theories to classical gravity solutions containing black holes. Thus, certain 
intractable strongly coupled field theory calculations are reduced to the study of black hole 
perturbations in General Relativity. 

One such example is the gravitational calculation of the ratio of the shear viscosity 77 to 
the entropy density s which yields the value ^ ^ | , and is generic for strongly coupled gauge 
theories with an Einstein gravity holographic dual (see e.g. Q]). On the experimental side, 
it has been argued that the value oirj/ s in the quark-gluon plasma produced in Relativistic 
Heavy Ion Collisions is indeed of the same order of magnitude (see e.g. jsj]). 

The ^ result can be obtained using a number of inter-related calculations. In jsl a holo- 
graphic calculation of the two-point correlator of the stress-energy tensor was performed and 
the shear viscosity was extracted from it using the Kubo formula. The transport coefficients 
can be extracted also from the linearized quasi-normal modes on the black hole background, 
which in the hydrodynamic limit, correspond to shear and sound modes of the gauge theory 

B. 

In j3] it was shown that the diffusive modes could be understood in terms of fluctuations of 



the black hole horizon using the language o 
a horizon with a fictitious viscous fiuid. In 



' the membrane parad 



igm js 



where one identifies 



9| this approach to calculations of shear viscosity 
and conductivity was formalized, showing how and when the horizon captures the response 
of the boundary theory. 

The hydrodynamics of non-conformal gauge theories is also characterized by a bulk viscos- 
ity coefficient C,. In general, holographic bulk viscosities do not exhibit the simple universal 
character of the shear viscosity. Also, the status of bulk viscosity in the membrane paradigm 
is more mysterious. In the old membrane paradigm Q], bulk viscosity is identified as neg- 
ative, with a ratio C,/ri = —1. Since the bulk viscosity of a real non-conformal fiuid must 
be positive (and generally also non- universal), it has been expected that a horizon based 
calculation will never correctly capture the bulk viscosity. The analysis of shear viscosity 
in {9] relies on the simplicity of the equation satisfied by the relevant transverse graviton 



mode on the black hole background. To compute the bulk viscosity, one must instead 
consider the trace mode h^. In the gravity solutions dual to non-conformal gauge theories, 
this mode may be mixed with fluctuations of scalar fields, and the calculations will be much 
more involved. 

The purpose of this paper is to study how the horizon captures the bulk viscosity in 
the framework of the fluid/gravity correspondence. In the conformal setup one considers a 
boosted black brane solution, characterized by uniform boost velocity u^, Hawking temper- 
ature T, and also, in general, chemical potentials and the corresponding charge densities 
p"-. This is dual to an equilibrium perfect fluid state of the boundary gauge theory. Next, 
one perturbs this solution by making the ansatz that the fluid velocity and thermodynamic 
variables are spacetime functions and solving the Einstein equations in a derivative expan- 
sion in Knudsen number i^fp/L <^ 1, where imfp is the mean free path and L the size 
of perturbations 10|. With the solution to first order in derivatives, one can compute the 
holographic boundary stress tensor at the first viscous order and read off the viscosities. The 
constraint equations projected on the boundary are the relativistic Navier-Stokes equations. 

In 111 it was shown that in these solutions the constraint Einstein equations projected 
on the horizon (Gauss-Codazzi equations) capture the same information in a simple way. 
In particular, the null focusing (Raychaudhuri) equation describing the evolution of the 
horizon entropy is equivalent to the viscous fluid entropy balance law. Using the zeroth 
order solution in derivatives one can read off the famous -j- ratio and thermal conductivities 



of holographic abelian and non-abelian gauge theories as well 12| . The various schemes are 
related by an RG flow in the radial direction [l^. 

In this paper we will extend this method to non-conformal fluids. We will show that the 
bulk viscosity arises from the flux of scalar fields across the null horizon and derive a novel 
and simple formula for the bulk viscosity in terms of derivatives of the horizon value (s) of 
the scalar field(s) with respect to thermodynamic variables such as the entropy and charge 
densities. 

We will consider ((i-|-l) -dimensional gravitational backgrounds holographically describing 
thermal states in strongly coupled d- dimensional field theories. The {d + l)-dimensional 
gravitational action in the Einstein frame reads 

1^ / V^d'-^'x (n-^ J2(^<P^f - v{<p.)^ + 1„ , (1) 
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where we set h = c = = 1 and canonically normalized the kinetic terms for the various 
scalar fields 0j. V{(f)i) represents the potential for the scalar fields and Igauge represents the 
action of gauge fields (abelian or non-abelian) A^. 

The null horizon focusing equation obtained by projecting the field equations of ([1]) on 
the horizon is equivalent via the fiuid/gravity correspondence to the entropy balance law of 
the fluid. Using this equation we will derive the following formula for the fluid bulk viscosity 

Here rj is the shear viscosity, s is the entropy density, are the charges associated with the 
gauge flelds AJ^, and (pf^ are the values of the scalar flelds on the horizon. We will apply 
([2]) to three general classes of gauge theory plasmas: non-conformal branes, perturbations 
of A/" = 4 supersymmetric Yang-Mills (SYM) and holographic models of QCD. 

The organization of this paper is as follows. In Section II we will briefly describe the 
horizon geometry and how the (Raychaudhuri) null focusing equation is mapped to the 
equation of the fluid entropy law. In a generic setup, we will show how the bulk viscosity 
at the horizon is encoded in the simple formula. In Section III we will apply the formula to 
the hydrodynamics of non-conformal theories such as those described by the near-horizon 
limit of Dp-branes, the Sakai-Sugimoto model in the quenched approximation, and non- 
conformal models that arise via a compatiflcation of a higher dimensional conformal theory. 



We then consider in 



The formula correctly reproduces the results of jl5| and jl6 |. 
Section IV deformations of A/" = 4 SYM. These include the M = 2* plasma studied in 
17 1, [l8|, a deformation of charged TV = 4 plasma studied in 19 1 and the softly broken 



conformal symmetry theories studied in 20| . Our formula will reproduce the results of these 
papers. In particular, we can obtain analytical results for the bulk to shear viscosity ratio 
without the need to resort to numerical methods. In Section V we will consider the bulk 
viscosity of non-critical backgrounds proposed as phenomenological holographic descriptions 



of QCD [2l|, |22[. We will show that our formula captures the bulk to shear viscosity ratio 
of these models in either the high temperature limit or in the adiabatic approximation, but 
otherwise deviates. In the last section we conclude with a discussion of the formula's range 
of applicability. 



II. BULK VISCOSITY FROM THE NULL FOCUSING EQUATION 



In this section we will start with a brief review of the horizon geometry and of the 



derivation of the null focusing (Raychaudhuri) equation ll|. Using the fluid/gravity corre- 
spondence, we map this equation into the fluid entropy law. We then show how the bulk 
viscosity is encoded in a novel and simple formula. 

A. Horizon geometry and the focusing equation 

We denote the coordinates of the bulk spacetime by = (r, x^), A = 0, ...,d. The x'^ are 
local coordinates on the horizon "H; r is a transverse coordinate, with r = on "H. Oat is a 
null covector tangent to the H. When raised with the bulk metric, it gives a vector field 

These choices fix the following components of the inverse bulk metric on "H: 

The pullback of qab into is the degenerate horizon metric 7^^. Its null directions are 
the generating light-rays of "H, i.e. ■y^^^'^ = 0. The Lie derivative of 7^1, along gives us 
the shear/expansion tensor, or "second fundamental form": 

Of,u = ^Ca^^ . (5) 
We can write a decomposition of into a shear tensor alu, and an expansion coefficient 6: 

d.u = + ^^7.. • (6) 
The expansion 6 can be expressed as 

e = v-^d.ivn = v-^d,S^ , (7) 

where u is a scalar density equal to the horizon area density, and S'^ = vl^ is the horizon 



area current 



It is convenient to raise indices with [G ^)^'^, which is the inverse of any matrix G^i, of 



the form 12 1 



G^. = Xl^u - b^K; b^i^ ^ . (8) 
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Here we introduced the scalar field A for later convenience: it will turn out that a matrix 
of the form ([8]) coincides at leading order with the metric h^v = rj^^ associated with the 
(typically fiat) metric of the hydrodynamic dual theory. 

Since 7^j/ is degenerate, one cannot use it to define an intrinsic connection on the null 
horizon, as could be done for spacelike or timelike hypersurfaces. The bulk spacetime's 
connection does induce a notion of parallel transport in "H, but only along its null generators. 
This structure is not fully captured by 7^,^; instead, it is encoded by the extrinsic curvature, 
or 'Weingarten map' O^*^, which is the horizon restriction of V^^^: 

e/ = v,t . (9) 



One can show that given an arbitrary G^y of the form ([8]), Q^j^ can be written as 12| : 



e/ = xe^piG-y^ + c/^ c/^ = K . (10) 

K is the "surface gravity", which measures the non-affinity of l^. The covector encodes 
the degrees of freedom in B^'' which are independent of 7^jy. In the hydrodynamics, these 
degrees of freedom will roughly correspond to the velocity and temperature fields. 
Re-expressing the identity, 

7^^iJ£^£^ = £^(VbV^ - VaV^j)^^ , (11) 

in terms of 6^*^ and imposing the Einstein field equation, we get the null focusing (Ray- 
chaudhuri) equation 

where Tab is the bulk matter stress tensor. 



B. Fluid entropy law and bulk viscosity 

In the following, we will consider gravitational solutions holographically dual to thermal 
states in strongly coupled field theories. We will begin with the action ([T]); its field equations 
read 

TZab - \9ABn =]^Y. {^^A4>i^B<P^ + \9ABV{<Pi)^ + Tf^'' . (13) 
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In (radially shifted) Eddington-Finkelstein coordinates the metric solutions that we are 
interested in take the form 

gAsdx^dx^ = dsfo) = -Crir + R)l^,lydx^dx'' + 2cR{r + R)l^dx^dr + c\{r + R)P^^dx''dx'' , 

(14) 

where at r = 0, the function ct vanishes and there is a horizon. The vector i^^ is a boost 
vector and has the dual role as a fluid velocity and the null normal to the horizon at r = 0. 
The indices /i, u etc. are raised with the fiat metric G^/^} = rj^y, and 

is the projection tensor. In addition there are the solutions to the scalar fields and gauge 
fields 

0,(r + i?), = A'^(r)5l, . (16) 

The Bekenstein-Hawking entropy density and the Hawking temperature associated with the 
metric are 

. = T=5r<^l;|M. (17) 

4 cr[R) 

From the solution for the gauge field, one can read off the values of the chemical potentials 
Ha and the charge densities p". The energy density and pressure associated with this solution 
are determined from the formulas 

de = Tds + /iarfp", e + P = sT + fiap"- ■ (18) 

With the thermodynamic quantities in hand, we can follow the ansatz of the fiuid-gravity 
correspondence, where one allows the thermodynamic variables, R, and i'^ to be functions 
of x'*. Now the zeroth order metric ( fT4l) and other fields as functions of x'^ are no longer 
solutions and must be corrected order by order in a derivative expansion. Let us see what 
we can learn from the Raychaudhuri equation (fT2l) . We write it in the following form 
corresponding to the action ([T]) 
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where we define the operator D = l^d^. Note, that does not appear in the equation 

since is null. In the equilibrium solution, the expansion and shear vanish identically, while 
the scalars and gauge fields parts vanish since these fields have only a radial dependence. 
Consider now the focusing equation at first order in derivatives of x^. In the case where 
T^y^^^ are general non-abelian gauge fields, one can show 
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12| that the first contribution 



appears at Oi^x^). Hence the focusing equation at order 0(x) reduces to 

^(1) = d^{st) = , (20) 

which matches the ideal entropy conservation law for the fiuid. Here the index on 9 refers 
to the order of the term in the derivative expansion. 

In proceeding to second order 0{x'^), one may worry that we might need information 
about the first order corrected metric. However, the structure of the Raychaudhuri equation 
is such that this is not the case. Consider first the shear term. The first order horizon shear 

-2 

can be calculated from the zeroth order metric. With the choice A = (4s) t*-! the shear 
squared term becomes 

where 

TT^, = p;p;9(,4) - -^p.^dp^' ■ (22) 

Next, the first three terms of f|T9|) to this order are: 

ra^^(^) + + + /c^o^^^^) _ ^^2^) . (23) 

However, at this order we can substitute the ideal equation 9^^^ = 0. Also, at first order 
there is always an ambiguity in the definition of our variables and s (or T). We remove 
this by the frame choice that always coincides with the entropy current, that is that 
5'(^)/^ = and s^^^ = 0. Therefore the 9^"^^ piece also vanishes. 

In 12| it was shown that any gauge field contribution at O^x"^) can be expressed in terms 



of (positive definite) entropy production terms associated with conductivities. From here on, 
we neglect this term and focus on the viscosities. Using k^^^ = 27rT, the focusing equation 
fITU]) has the form 

u^n = \d,s^ = ^^^'^^"^ + ^ +■■■ (24) 
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where the dots indicate the gauge field terms. Notice that the scalar term is a derivative of 
0j-''''(r, x'^). We take the derivative of this solution with respect to x'^ and then evaluate at r = 
0. This is equivalent to D(j)f{x^) of the horizon values. The horizon values (j)f are functions 
of the entropy density s{x^) and charge density (or the Hawking temperature T(x^) 

and chemical potential /ia(x'^)). Thus we can write 

In the last equality we used the ideal conservation law (1201) and also the ideal conservation 
equation for any charge current present 

d^J^"" = d^{p''n=^ ■ (26) 
Therefore, the Raychaudhuri equation to second order takes in the fiuid/gravity correspon- 



dence the form of a fiuid entropy balance law 2^ 



5.(^n = fw''^ + ^(5pn', (27) 

with the standard shear viscosity obeying rj/s = l/47r, but now with an additional bulk 
viscosity term. Substituting f l25|l into flM|) we find the bulk to shear viscosity ratio implied 
by the Raychaudhuri equation is simply 

or equivalently 



C 
s 



In the absence of charges we can use the speed of sound = and recast the 



formula in the form 



r^-'-i^y- (30) 



III. NON-CONFORMAL BRANES 



In this section we will use the formula (|2]) to calculate the bulk viscosity of non-conformal 
hydrodynamic models realized holographically as gravitational backgrounds with scalars, 

9 



dua' 

Q 



to non-conformal branes. We will consider three examples: Dp-brane here p 7^ 3 



15| . the Sakai-Sugimoto model in the quenched approximation 15| and non-conformal 
brane backgrounds obtained from the dimensional reduction of pure gravity 16 1. 

We will show that in all the three models, the horizon formula reproduces correctly 
the previous results obtained by either the holographic hydrodynamics or by a holographic 
calculation of the stress-energy two-point correlators and the use of the Kubo formula (or 
the equivalent quasi- normal modes analysis). 

A. Dp-branes 

The decoupling limit of stacks of coincident Dp branes with p 7^ 3 yields a gauge/gravity 



correspondence of non-conformal gauge theories [25|. The (super)gravity solution is confor- 
mal to AdSp^2 x S^~^, and a key feature of the solution is the existence of a non-trivial 
scalar dilaton field. 

The {p + 2)-dimensional action for the metric and dilaton, reduced from the ten- 
dimensional supergravity action reads 

s . J- /,.«.^ (n - \ie,f - (Iz£Ki^,^*") , (31) 

where we canonically normalized the scalar kinetic term. The zeroth order equilibrium 
solution for the metric components can be put in the form of the metric (!T4|) with 

2 \ 2 / \ \^ (p-3)(p-6) 

Crp{r) = r p f{r); Cx{r) = [r + R) p ; cj^(r) = r , (32) 

where /(r) = (1 — ;;:^)^~^. For the scalar dilaton, the solution is 

0(r) = ^{^-p) (^] ' ln(r + R) . (33) 



2 \ P 

The thermodynamics of this equilibrium solution is characterized by the Hawking tem- 
perature and the Bekenstein-Hawking area entropy density 

T = — = l^Ri^~p)/^- s = - = (l/4)i?(9-P)/2 , (34) 

The value of the dilaton on the horizon as a function of the entropy density is 

_ 1 
<P^ = ^{3-p) (^^) ' ln((4.)^) . (35) 
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Taking the derivative of cj)^ with respect to the entropy density, we find using ([2]) 

e 2(3 



(36) 



?7 p(9 — p) 

which agrees with jl^, who computed the ratio using the hnearized quasi- normal modes 
analysis of this background solution. 



B. The Sakai-Sugimoto model 



26| proposed as a 



We now consider the hydrodynamics of the Sakai-Sugimoto model 
holographic model of QCD with fundamental fiavours. We will consider the quenched ap 
proximation discussed in 15|. The effective five- dimensional action reads 

1 f . — /_ 1_„., 1 



S 



-gd^x [n-\{dff - ]^{d^f - ^{dw) 



V 



(37) 



167r i " " V 2 
where P is a potential function that depends on the three scalars f,^,w. The zeroth order 



solution again can be put in the same form as (fT4|) . The entropy density is 



1 ^5 



(38) 



and the zeroth order scalar field profiles are 



$ = - ln(r + R), 



w 
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4^15 



ln(r + R) 



(39) 



Evaluating at the horizon and using the formula we find 



^ _ 4 / 9 10 169 
^~25Vl6^25^240 



4 
15 



(40) 



which agrees with the result found in 15[ |. 



C. Non-conformal branes 

It has been demonstrated in |l^ that solutions which asymptote locally to non-conformal 
branes can be obtained from higher- dimensional asymptotically AdS solutions via a dimen- 
sional reduction and a continuation in dimension. This generates a class of holographic 
descriptions of rf-dimensional non-conformal hydrodynamic models parametrized by posi- 
tive number a, where the energy density e and the pressure p are related by the equation 

11 



of state e = (2cr — where 2(t ^ d (recall that e = {d — l)p is the equation of state for 
conformal hydrodynamics). By using the hydrodynamics expansion, one can find a general 



formula for the ratio of the bulk to the shear viscosity in this class of models, it reads 16 1 

C 2{2a-d) 



r] (d- 1)(2(7-1) ■ 
We will now derive this result using 

In 161] the authors work in the dual frame, where the action is 



(41) 



S=^l d'+'x^ge^ [u + '^\J_/ id<pf + 2a{2a - 1)) . (42) 

The corresponding black brane solution reads 

ds^ = (r + RfP^^dx^dx'' - firy^tdx^dx" + 2i^dx^dr (43) 
(f)={2a- d) ln(r + R) , (44) 

where /(r) = 1 - (^)"^" and R=^- 

The entropy density in the dual frame is not the area density v, instead we have 

s = ^exp(0^)w, (45) 

which for this solution is 

s = Ir'"-' . (46) 

Our analysis of the Raychaudhuri equation in Section II B depends on the entropy density 
being proportional to an area density. Therefore, to apply our formula (Ej), we need to 
make a conformal transformation of (jH]) to the Einstein frame. The bulk viscosity of course 
does not depend on the choice of conformal frame. In the Appendix we re-analyze the 
Raychaudhuri equation in the dual frame and arrive at the same result. 

To transform the action fH2|) to the form of ([1]), we first make the conformal transforma- 
tion 

20 

9ab = ed-^gAB, (47) 

where the tilde represents the metric in the Einstein frame, and then canonically normalize 
the scalar field. On the horizon, the scalar field solution then has the form 



^» = ^hl<2^fM^HIi). (48) 

y/d — 1 

Re-expressing this in terms of the entropy density, (l46|l . and using the formula (j2]) we obtain 



(j^Tj) in agreement with the result found in 16 1. 
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IV. PERTURBATIONS OF 7^ = 4 PLASMA 



In this section we consider deformations of A/" = 4 SYM. These include the M = 2* 
plasma studied in [l^], [l8|, a deformation of charged A/" = 4 plasma studied in [l^ and 
the softly broken conformal symmetry theories studied in [2^. We will show that the bulk 
viscosity analytical formula agrees with the previous calculations, some of which required 
previously a numerical analysis. 



A. J\f = 2* plasma 

The M = 2* four- dimensional gauge theory is defined as a mass deformation A/" = 4 
SYM, where the bosons and fermions masses and mj break the conformal symmetry. 
The M = 2* gauge theory has, in the limit of large t'Hooft coupling, an explicit supergravity 



dual description called the Pilch- Warner flow 



271 ]. The five- dimensional geometry dual to a 



finite temperature state is characterized by three independent parameters: the temperature 
and the horizon values of two scalar fields a and x- 

The effective five- dimensional action for this model is given by [Ij 

S = ^l V^d'x (n - ^{daf - ^{dxf - 4P^ , (49) 
where we canonically normalized the scalar fields a and y and V is their potential. The 



values of the scalar fields at the horizon are denoted by [171] 



a = Si] x = ^2] cx = h ■ (50) 



In the high temperature regime of the M = 2* black brane, the 



gravity field equations can 



with 



be solved analytically in an expansion in parameters (5^ <^ 1 |17 . 

1 /m,x2 v/2r(3/4)^ m/ . n ^ ^ 2 ^ A^^ 



'2Atx 

The entropy density associated with this solution is s = 5|/4. Plugging the derivatives of 
the horizon values of the scalars with respect to the temperature in ( 130|) we get 

^ = vt{72x + Qy) , (52) 



where we defined 

_ r 
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Since we work at linear order, we can use the zeroth order (conformal) value for the speed 
of sound, V. 



and arrive at 



c 

7] 



Comparing to the parametrization C/v = n ^ ^Pf V i^i Um, we see that /3j 



This is consistent with the values reported in 18|, where f3^ ~ 8.001 and f3^ ~ 0.6666 were 



(54) 



8 and l3j = 2/3. 



obtained via numerical analysis of quasi-normal sound modes. 



B. Deformation of charged J\f = 4 plasma 

The previous examples that we considered did not involve charges. An interesting test 
case for our formula is the mass deformation of the charged A/" = 4 plasma studied recently 
in 



19] . The action for this theory is given by 



s = — [ d^'x^ (n- -0^/3^2 - -0-2(90)2 + 402/3 + 80-1/3 + sl] 

167r y \ 4 3 / 



(55) 



where 6L represents the conformal symmetry breaking mass deformation due to a second 
scalar field x, 

SL = -I{dxr-'^X' + --- (56) 

The idea of the calculation is to work perturbatively in the x field. One makes the ansatz 
that X = -^Xi) where A is the coefficient of the non-normalizable mode of x- The gauge field 
Ab and metric functions receive corrections at 0(A2). At zeroth order in A (the unperturbed 
solution), the thermodynamic variables are defined in terms of constants k and /3: 

^ K + 2 /3 2txT ^2 

Jk{1 + k) Vl + K /tt ^JK 



The ratio of the temperature to the chemical potential obtains a minimum at k = 2, where 
the black hole undergoes a second order phase transition. The scalar field has the horizon 
value 



K . (58) 



The zeroth order solution should have a vanishing bulk viscosity, but the scalar field has 
a non-trivial horizon value, which seems to imply a potential contradiction. Let's consider 
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the general formula ([2]). The Eq. (2.47) in 19| for the entropy and charge density imply 
that 

^2 



K 



(59) 



47r2 

Using fl59p we can express 0^ as a function of p and s, and it is straightforward to see that 
the formula for the bulk viscosity ([2]) yields C = as expected from conformal invariance at 
this order. 

Non-zero bulk viscosity arises at the order O(A^) from the horizon value of the x held, 
where the parameter A is related to the fermionic mass M and the temperature by 

,_(«; + 2) M 

23/4(1 + ^)3/% r" • ^ > 

Thus the expansion in A is again an expansion in the limit of high temperature, where 
T ^ M. The field does not contribute at this order since its form changes at order O(A^) 
and therefore its contribution to the bulk viscosity is at order O(A^). 

We can compare to the numerical calculation done in 19| at the critical point (k = 2). 
The field x takes the form 



X 



,H 



Ac! 



h(Q) 



(61) 



The coefficient cj^°^ = 0.7464562054847809 (see Table I in [l9|]) is a pure number obtained 
numerically by solving ordinary differential equations subject to boundary conditions at the 
horizon and at infinity. 

When using the formula (|2]) we hold k fixed and vary the temperature T with respect to 
the entropy density s and the charge density p. Using the relation between T, s, p in Eq. 
(2.47) in [l9| is easy to see that ([2]) gives 



c 

7] 



cT^fX^ ~ 0.0619A^ 
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(62) 



which is consistent with the numerical result of ~ 0.06218A^ obtained by analysis of the 



quasi- normal sound modes 19 1 



C. Softly broken conformal symmetry 



A general formula for bulk viscosity of theories, where the conformal symmetry is softly 



broken has been derived in 



20| . One considers a bulk gravitational action with a scalar field 



15 



(63) 



The conformal symmetry of the field theory is broken by a relevant operator of conformal 
dimension A. In the duality prescription, the field is dual to Oa', the expectation value 
of Oa and the source term A'^"^ can be read off from a boundary series expansion of </>. 

In the high temperature A/T ^ 1 expansion, the relevant physics is simply captured by 
the scalar field propagating on the fixed black brane background. For our purposes, all we 
need is the (zeroth order) scalar field solution 

0(r6) = (A6)'^^^-;^P-.i+a(-1 + 2(r6)'^) , (64) 

where P is a Legendre function of the first kind and b = d/{47TT). 



To compare more easily with the results in 20|, we will work with the bulk viscosity 
directly. Taking the derivative of with respect to the temperature, evaluating at the 
horizon, r = b~^, and using our formula in terms of the temperature derivative ( l30l) . we find 

The entropy density and speed of sound for background are simply (in units where = 

(167r)-i) 

Using these results and the F function identity, 

T{z)T{z + 1/2) = 2i-2^v^F(2z) , (67) 

we find 



AtiTJ ^ (47r)2('^-^) (d- 1)2 \^F(A/rf+ 1/2) / \T 

which matches Eq. 4 of 201 ]. 
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V. HOLOGRAPHIC QCD MODELS 



In this section we will consider holographic phenomenological models of QCD, based on 
gravitational backgrounds with a scalar. In this setup one starts with one scalar field coupled 
to gravity in five dimensions with the action 



One considers the scalar field $ = </'(r), with the metric ansatz 



(69) 



(70) 



29| one makes the gauge choice that the scalar field itself is the radial coordinate, 



Following 
i.e. $ = r. 

With the choice of metric ansatz, one finds the entropy density 

g3A($H) 



(71) 



which implies that 



and the formula ([2]) gives 



dcf)" _ 1 



(72) 



r] 9 

We can go a step further by using the rr metric and the scalar field equation 



(73) 



6Ah' + h{2AA" - 1) + 2e'''V = 0, 

AA' -B' + - + —V = 
h h 



(74) 
(75) 



evaluated at the horizon {h = 0). These imply that A'{^h) = V{^h)/'^V'{^h), and thus 



y($ 



(76) 



This result agrees with the ratio found in 



21 



22| using the holographic Kubo formula 



rj V^($h) 



(77) 
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only when the overall factor Cb = 1. To explain the origin of this factor, we briefly review 



the analysis of 



21 



22| . which starts by considering the Kubo formula 

C = -^lim-/mG^(w,0) , (78) 

9 t^-S'O oj 



where is the Fourier transform of the retarded Green's function of the trace of the gauge 
theory's stress tensor 

G",., k) ^ -. I #.*e-'--.-.«) < m, T/(0. 0)1 > . (79) 

In the AdS/CFT prescription, two point functions of stress tensor operator T^i, can be 
calculated using the metric perturbations -f^^v Using 50(3) invariance and the gauge $ = r, 
one can show that the relevant metric perturbation dual to Tj* is 

6g = diag{gu, gn, gii,gii, gw) , (80) 

where, to linear order, 

2B 

gu = -he'^^il + Hu); gn = e^^{l + H^,); g^ = —(1 + H„) . (81) 

Using the ansatz of harmonic time dependence for the perturbations, e~*'^*iJ^jy(<l'), one can 
show the equation for Hu decouples 

To find the bulk viscosity it is sufficient to solve for Hn. This is done subject to the 
boundary conditions that Hu — > 1 at the boundary and that it is purely infalling at the 
horizon ($ — > $//): 

Hn^c,{^H-^)~^, (83) 

where is some normalization factor. With the solution in hand, one can show that the 
Kubo formula reduces to just fl77|) . 

Determining the bulk viscosity via the Kubo formula in this class of models amounts to 
determining by solving numerically fl52]) with a; = and the boundary conditions imposed 



at the horizon and at infinity. As was shown explicitly in [2l|, |22|, ci, is not 1 in general and 
therefore there is a discrepancy between the Kubo formula answer for the bulk viscosity and 
the value computed with the focusing equation. A plot of Cb as a function of \h = e^" is 

18 



presented in Figure 5 of 22|], which shows that the focusing equation agrees with the Kubo 
formula only in certain limits. One feature of this plot is that Ch approaches unity in the 
high temperature regime. In this case the horizon is essentially near the boundary. 

Another important approximation where the horizon does capture the bulk viscosity is 
in the adiabatic regime discussed in 29|. Note that when the potential is an exponential 
function of the scalar field (see, for example the Dp branes in Section III) the speed of sound 
is a constant related to V'{^h)/V{^h)- In the adiabatic regime, one allows V /V and the 
speed of sound to be a slowly varying functions of the horizon value of the scalar field. In 
terms of the entropy density and temperature, this translates into 

lns = - d^— H 

InT^J d*(-----)+... (84) 

where the ■ ■ ■ represent the slowly varying terms in In 
the fluctuation equation ( 1821) for Hn simplifies in adiabatic regime. As a result, the solution 
for Hii is independent of $ and the boundary conditions require Cb = 1. 

The trivial radial behavior of the u = linearized fluctuation equation in the adiabatic 
regime is reminiscent of the behavior of the equation discovered by Iqbal and Liu 9j that 
governs the radial flow of a shear viscosity defined on timelike surfaces of constant r between 
the horizon and boundary. In the low frequency limit, the radial dependence in this equation 
always drops out and it is convenient to evaluate the shear viscosity of the boundary gauge 
theory on the horizon. Our results suggest that a similar equation for the bulk viscosity 
would only allow one to work at the horizon if the additional condition that the potential 
and thermodynamic variables are slowly varying functions is imposed. 



22| it was demonstrated that 



VI. DISCUSSION 



Using the null focusing equation, we derived a simple novel holographic formula ([2]) for 
the bulk viscosity of field theory plasmas at strong coupling. The formula was expressed in 
terms of the dependence of scalar fields at the horizon on thermodynamic variables such as 
the entropy and charge densities. We applied the formula to three classes of gauge theory 
plasmas: non-conformal branes, perturbations of A/" = 4 SYM and holographic models of 
QCD. In the first two classes of models the formula agreed with previous calculations based 
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on quasi-normal modes analysis (or equivalently the holographic Kubo formula) and the 
hydrodynamic expansion. Note, however, that in all these cases, the calculations of the 
bulk viscosity in the literature where performed essentially in the high temperature approx- 
imation. In the last class of phenomenological holographic models of QCD, we found an 
agreement in two approximations: the high temperature limit and the adiabatic approxima- 
tion. 

Let us discuss the range of applicability of the formula. The null focusing equation 
encodes dynamics of the horizon. For the boundary field theory this is the IR regime. The 
bulk viscosity has in general a non-trivial dependence on the energy scale. For instance, 
in QCD type theories, dimensional transmutation expresses this non-trivial dependence by 
relating the running coupling g at the temperature scale T to the dynamical scale A 



T 



e ^ , (85) 



where h is the one- loop /3-function coefficient. In such theories, the IR horizon calculation 
is not expected to capture this running. 

In the high temperature limit, the horizon is basically near the boundary and essentially 
captures also the UV physics; this is the reason for the agreement between the different 
types of calculations of the bulk viscosity. In the adiabatic limit on the other hand, when 
considering the low frequency limit, the radial dependence in the calculation of the bulk 
viscosity drops out and horizon and boundary calculations yield the same result. It would 
be interesting to work in the membrane paradigm formalism to see if this is indeed the case. 



A bound on the bulk to shear viscosity ratio has been proposed in 30| : 



A y2 

d-l 



> 2 . (86) 



Interestingly, the bound is satisfied by the formula ([2]) in all the cases that we studied in the 
paper. In particular, the bound is satisfied in the high temperature expansion used to study 
the various perturbations of the TV = 4 SYM, while in the adiabatic regime, the formula ([2]) 
saturates the bound. This follows from (!77|) with Cb = 1, together with the approximation 
for f ^ 



s ■ 
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APPENDIX 

Here we show explicitly that the choice of conformal frame in Section III C does not effect 
the bulk viscosity. The equations of motion for the dual frame action fH21) imply that 

R,ui^t = -^{D^f + , (87) 

Za — a 

so that the Raychaudhuri equation has the form 

-Fd^e + t,e- _ ^ _J_(^jj^)'2 + (88) 

a — i la — a 

Note that 9 is expressed in terms of the area current as before, v~^d^j_{vl^). Here, instead 
we want to work with the entropy current = (l/4)e'^f£'^. One can define a modified 
expansion 

e = {e'^v)-^d^{e^vF) =D(j) + e . (89) 

The entropy current S^s'^ = s6. Now let's consider flSSl) at lowest order. First, we rewrite 
the 2nd covariant derivative of (p in the following way: 

rrv^v^0 = rv^(rv^0) - KD(f). (90) 

The first term is now O(e^), but the second is first order, by virtue of the geodesic equation 
ii'V/'' = Kt. Thus we are left at 0(e) with 

k{9 + Defy = Ke = 0. (91) 

So, as expected at ideal order, the Raychaudhuri equation is equivalent to the vanishing of 
the entropy current d^{si^). 

Now let's consider the equation to O(e^). We write it in the following form 

27rT9 = D9 + -^6^ + 7r^,7r^^ + -^—-{D<pf + D{D{<p)) . (92) 
a — i 2a — a 
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Now one can impose the ideal equation 9 = —Dcj) in the O(e^) terms on the right hand side. 
Expressing the 1)0 term again in terms of derivatives of the thermodynamic variables, we 
find 

Using the scalar field solution, 

(f)^ = {2(x - d) ln(4s)^ , (94) 
one can read off a bulk to shear ratio that agrees with (14T|) . 
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